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1.  INTRODUCTION 

Many  treatments  of  the  acousto-optic  effect  — 
the  interaction  of  light  with  sound  —  have  been 
published.' 6  7 *  These  treatments  assume  an  infinite 
slab  of  isotropic  material  in  which  an  acoustic  wave 
is  propagating  parallel  to  the  slab  surfaces.  A 
coherent,  monochromatic  electromagnetic  plane 
wave  is  incident  on  the  medium  from  one  side,  and 
emerges  from  the  other  side  broken  up  into  discrete 
diffracted  plane  waves  of  all  integral  orders  (fig.  1). 
The  acoustic  disturbance  modulates  the  refractive 
index  of  the  medium,  producing  what  may  usefully 
be  visualized  as  a  diffraction  grating. 

In  previous  work,  the  researchers  have  used  a 
scalar  wave  theory  (neglecting  polarization  effects) 
and  have  (except  for  Berry4)  assumed  that  re¬ 
flection  from  the  two  interfaces  can  be  neglected. 
Berry  has  considered  reflections  (in  a  scalar 
theory)  and  shows  that  reflected  diffracted  beam 
intensities  for  a  slab  of  a  certain  thickness  are 
directly  proportional  to  the  transmitted  diffracted 
beam  intensities  for  a  slab  of  twice  that  thickness. 
(Actually  he,  like  the  other  authors,  had  in  mind  a 
liquid  medium.)  The  proportionality  constant  is 


For  liquids,  this  quantity  tends  to  be  small  (0.01, 
or  so),  whereas  for  solids  it  may  be  appreciably 
larger. 
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In  an  attempt  to  model  the  acousto-optic  effect 
in  a  transparent  solid,  the  author  has  undertaken  to 
generalize  previous  results  to  include  reflected 
waves,  and  also  to  consider  anisotropic  media.  The 
work  reported  here  uses  a  scalar  theory  for  an 
isotropic  medium  (like  earlier  work),  but  with 
rigorous  boundary  conditions  —  i.e.,  solving  for 
reflection  as  well  as  transmission.  The  problem  is 
treated  by  using  a  series  expansion  in  the  parameter 


(where  S  is  the  intensity  of  acoustic  modulation, 
and  A  and  X  are  acoustic  and  optical  wavelengths 
respectively);  hence,  the  expansion  converges  only 
when  this  parameter  is  small.  A  vector  treatment  of 
the  interaction  in  an  anisotropic  medium  will  be 
reported  in  a  subsequent  paper. 

Motivation  for  this  research  is  given  by  the 
experimental  work  of  Berg,  Lee,  and  Udelson,*  9 
who  are  investigating  the  acousto-photorefractive 
effect  in  y-cut  slabs  of  lithium  niobate,  an  aniso¬ 
tropic  material.  In  their  work,  the  acoustic  dis¬ 
turbance  is  given  by  a  surface  acoustic  wave 
(SAW)  excited  by  interdigital  transducers  and 


ACOUSTO- OPTICAL  INTERACTION 


Figure  I.  Incident  wave  is  split  by  acoustic  wave 
into  discrete  diffraclrd  orders. 
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propagating  along  the  c-axis.  The  slab  is  subjected 
to  a  high-intensity  short-duration  laser  pulse, 
which  produces  a  “freezing”  of  the  instantaneous 
refractive  index,  an  effect  which  can  be  used  in  a 
memory  correlator.  Also,  work  is  being  done  on 
real-time  acousto-optic  correlators  and  convol¬ 
vers  using  two  SAW's  simultaneously.  For  a  com¬ 
plete  understanding  of  these  processes,  a  thorough 
theoretical  treatment  of  the  diffraction  is  needed. 
This  paper  treats  the  isotropic  time-independent 
case,  and  can  be  regarded  as  a  preparation  for 
treatment  of  the  full  time-dependent  anisotropic 
case.  (The  time  dependence  is  a  trivial  addition 
giving  a  frequency  shift  of  the  diffracted  waves. 
Anisotropy,  however,  introduces  considerable 
complication,  because  of  the  acoustic  coupling  of 
ordinary  and  extraordinary  modes.)  This  paper 
does  not  treat  the  spatial  dependence  normal  to 
the  surface  (exponential  falloff). 


2.  THEORY 

The  physical  system  to  be  analyzed  is  an 
infinite,  isotropic,  and  lossless  slab  of  thickness  L. 
An  acoustic  disturbance  is  propagated  through  the 
medium  in  the  z-direction  parallel  to  the  slab 
surfaces,  and  the  x-direetion  is  chosen  normal  to 
the  slab  surfaces.  A  plane  wave  of  monochromatic 
light  is  incident  in  the  x-z  plane,  at  an  arbitrary 
angle  and  linearly  polarized  in  the  y-direction. 

The  motion  of  the  acoustic  disturbance  is 
neglected  and  its  effect  is  incorporated  as  a  varia¬ 
tion  in  the  dielectric  property  of  the  material 
according  to 

e  —  e  +  S  cos  qz  , 
where 

f  =  dielectric  constant  of  unmodulated 
material. 

S  =  modulation  amplitude,  and 


2rr 


This  assumption  corresponds  approximately  to 
the  interaction  of  a  pulsed  light  wave  whose 
duration  is  short  compared  with  the  acoustic  wave 
period.  The  wave  equation  for  a  light  beam  of 
radian  frequency  at  and  amplitude  E  is  then 

—  ^V:  H — fe'j  E  =  ficosqz^E  .(1} 

where  c  is  the  velocity  of  light. 

A  propagating  electromagnetic  mode  inside  the 
medium  can  be  resolved  into  plane  waves  as 
follows: 

E  =  £  En6  exp(ikx)  exp(iy„z)  As  .  (2) 
ns 

where 


y  =  z-component 

of  incident  light  wave 

vector. 

Yn  =  y+nq 

(3) 

“< 

8  Wn 

II 

r* 

IK 

2 

.  CD 

A  =  S  —  = 

c  q 

4  v 

Note  that  this  is  not  a  power  series  expansion  in 

A,  since  k  =  k(A).  and  therefore  it  is  not  unique. 
It  can  more  aptly  be  termed  a  plane-wave  expan¬ 
sion.  Substitution  of  equations  (2)  and  (  l)  into  the 
wave  equation  ( 1 )  yields 

£  +  Yn~  +  ZM:n,|A'  j 


exp(ikx)exp(iynz)  A*  = 


exp(ikx)exp(iy„z)A“  (<>l 
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The  zero-order  component  E^  represents  a 
"normal”  mode  of  propagation  for  the  medium  — 
i.e.,  a  solution  of  the  homogeneous  equation 


From  this  we  see  that 


1 


c 


2 

e-y 


Equating  like  powers  of  A  in  equation  (6)  yields  the 
recursion  relation  in  Ens 


n{n  4-  2a)Eng  =  '/j  (En-i.„-i  +  En+i.*-!)  - 


where 


— r  2  KlEn.»-l 

q  1=1 


(7) 


Given  the  angle  of  incidence  of  the  light  wave,  ft  y 
is  determined  by  y  =  to/c  sin  ft  and  equation  (7) 
generates  the  corresponding  propagating  optical 
mode  of  the  perturbed  medium. 


The  above  solution  can  be  written  as 
E  =  exp(  ikx)e(z),  where  e(z)  will  be  shown  to  be  a 
Mathieu  function.  Equation  (1)  becomes 

f-r-y  +  -y  t  -  k'  +  -7  J  cos  zq]  e(z)  =  0 
|_dz  c  c  J 

This  is  Mathieu's  equation: 

J*f 

— 7  +  (a  -  2h  cos  2{)f  =  0  , 

d< 

where 


kJj  =  (2cr)*  +  O(A')  . 

-2A 


2<  =  q* 


c  q 


Thus 

e(z)  =  £  Ens  exp(iynz)As  =  f(a,<)  , 

ns 

where  f(a,£)  is  a  Mathieu  function.  This  function 
will  not  in  general  be  one  of  the  well-known 
tabulated  periodic  functions  ce„  and  se„,  but,  in 
general,  will  be  a  Floquet  solution, 10  of  the  form 
eM*  0(z),  where  0(z)  is  periodic.  However,  when 
a  =  N/2  (N  an  integer,  the  Bragg  condition),  the 
solution  is  a  linear  combination  of  the  cen  and  se„ 
(see  appendix  A). 

The  preceding  paragraphs  characterize  the 
natural  propagating  modes  for  an  infinite  sinusoid¬ 
ally  modulated  medium.  For  a  Finite  thickness 
slab,  0  <  x  <  L,  the  electromagnetic  boundary  con¬ 
ditions  require  that  the  solution  be  an  expansion  in 
a  subset  of  these  natural  modes,  requiring  a  new 
index: 


E  = 

2  Emn,  exp(ikmx)exp(iynz)A,‘ 

mns 

(8) 

Also  a 
mitted 

backward-travelling,  reflected,  and  a  trans¬ 
wave  are  required: 

F  = 

2  Fmn*  exp(— ikmx)exp(iynz)As 

mns 

,  0) 

R  = 

2  Rn»exp(->Pnx)»**p(iYnz)As  , 
ns 

(10) 

T  = 

2  Tns  exp(ipnx)exp(iynz)As  , 

ns 

(II) 

where 

p£ 

__  &  2 

2  Yn  • 

(12) 

c 


Now  the  wave  equation  ( I )  yields  a  new  recursion 
relation: 

(n  -  m)(n  4-  m  4-  2cr)Emn„  = 


m.n~  l,s~  1 


4-  E 


m.n+ 1.*- 1 


)~ 


V  W  MrLtrhian,  7Tteor\  and  -t/tpfirdNoft  uf  MntHtru  functions.  < >tfnrd 

/Wmm, 
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2  X  KmtEm,n,s-i  (13) 

q  X  =  1 


There  is  a  similar  recursion  relation  for  F. 

The  boundary  conditions,  applied  at  x  =  0 
and  x  =  L,  guarantee  continuity  of  tangential  E 
and  H: 


Enns  "F  Fnns 

—  £  (^mns  d"  Emns)  i  (Ida) 

m 

(»*n) 


Pn^ns  d"  knuEnns  knl|Fnns 

s 

,  Z  (km.s- 1 Emn  t  -  kmns-  X  Fmn  f  )  ,  (14b) 

m  1=  0 
(*ns) 


Pn«Enns  d  Mm,F nns  exp(ipnL)Tna 


—  £  51  (Pm.s-lEmnl  d- Mn  s_|  Fmn £  )  ,(14c) 
m  i=  0 
(s^ns) 


kn<> EniiEnn»  kn0  Mn0  Fnns  pn  exp(ipnL)Tns 

S  9—  X  —  1 

-  £  Z  I 

m  X  =  0  j  =  0 

<*n) 

m.s— X  Pm,  X  -  jEmnj  kmg-(Mm,t— jFmnj)  , 

(14d) 

where 

exp(ikmL)  =  £  PmlAl 
X  *  0 


exp(-ikmL)  =  £  Mm£Al  (15) 

x  =  o 

If  the  incident  plane  wave  strikes  the  x  =  0 
interface  with  a  z-direction  cosine  of  y„  then  the 
aero-order  terms  Eom,  Fooo>  poo-  and  TM  are  just 


those  given  by  the  familiar  Fresnel  relations  for  an 
unperturbed  medium.  The  first-order  terms  are 
then  generated  from  the  recursion  relations  (which 
yield  E0,„  E0_„,  F0i„  and  F0_„)  and  the  boundary 
conditions  (which  yield  R,,, E,,„  F,,,,  T,,  and 
R_ ,  „  E_ ,  _  F_ ,  _  1 1,  T_ , ,)  This  process  can  be  con¬ 
tinued  indefinitely,  and  is  an  unambiguous  way  of 
generating  the  solution  for  a  given  incident  plane 
wave,  provided  A  is  “small,”  so  that  the  series 
converges.  The  recursion  relation  ( 1 3)  serves  not 
only  to  generate  the  Emns,  but  also  to  generate  the 
Km{.  When  m  =  n,  the  left-hand  side  of  (13)  is 
clearly  zero.  The  Km-S_m  is  determined  in  terms  of 
low-order  values  of  Kmt  and  Emn8.  (Enns  and  Fnns 
are  determined,  as  previously  stated,  by  the 
boundary  conditions.) 

The  left-hand  side  of  equation  (13)  is  also  zero 
whenever  n  +  m  4-  2a  —  0.  This  is  a  funda¬ 
mental  problem  and  in  fact  invalidates  the  pro¬ 
cedure,  but  it  occurs  only  when  a  =  N/2  for 
some  integer  N  —  the  Bragg  condition.  This 
difficulty  can  be  dealt  with  in  practice  by  simply 
treating  an  approximate  problem  a  =  N/2  +  £. 
An  elaboration  of  the  theory  is  given  in  appendix  A. 

3.  GEOMETRICAL  INTERPRETATION 

There  is  a  simple  geometrical  construction 
which  is  a  useful  aid  in  visualizing  the  above 
treatment  of  the  acousto-optical  interaction.  This 
construction  derives  from  one  often  used  to  de¬ 
monstrate  Snell's  law.  Thus,  in  figure  2,  two  views 
of  the  simple  unmodulated  boundary-value  prob¬ 
lem  are  given:  (a)  is  self-explanatory;  (b)  is  more 
abstract.  Here  the  circles  (in  wave-vector  space) 
represent  the  naturally  propagating  modes  in  the 
various  media.  Their  diameters  are 

to  to  -  -  to 

—  ,  —  V  f  -  aft*  — 

c  c  v  c 

Figure  3  represents  the  acoustically  modulated 
medium.  Here  are  represented  the  higher-order 
(m,n)  terms  —  plane  waves  generated  by  the 
recursion  relation  (13)  and  the  boundary  condi¬ 
tions  (14).  The  (0.0)  vectors  are  as  in  figure  2. 
Equation  (13)  generates  (0.1)  and  (0.— 1).  Then 
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(a)  PHYSICAL  SPACE  (10  WAVE-VECTOR  SPACE 


Figure  2.  Representation  of  zero-order  interac¬ 
tions  (Snell's  law)  in  physical  and  wave-vector 
space;  boundary  eonditons  require  that  iu\  c  c  sin  0 
be  constant  for  all  media. 


PLANE-WAVE  DECOMPOSITION 


Figure  3.  Internal  backward  and  forward  travel¬ 
ling  waves;  to  each  “lattice  site”  in  w are-rector 
space  corresponds  a  complex  amplitude,  whose 
value  is  given  by  the  recursion  relation  (eq  (13)) 
and  by  the  electromagnetic  boundary  condition. 

equation  (11)  generates  (1.1)  and  ( —  1 . —  1 ).  Then 
equation  (13)  generates  (1.2)  (1.0).  (0.2).  (0.— 2). 
(—1.0).  and  (  —  1.— 2).  Equation  (14)  generates 
(2.2).  (—2.-2),  and  a  second-order  term  for  (0.0). 
and  so  on.  It  will  he  recognized  that  the  recursion 
relation  imposes  a  vertical  coupling  on  the  diagram, 
while  the  boundary  conditions  impose  a  horizontal 
coupling. 

It  is  important  to  note  that,  in  figure  3.  the 
"circle"  is  not  a  true  circle,  hut  is  deformeii  by  the 
A-dependence  of  the  km. 

The  Bragg  problem  (o  =  N/2)  can  be  seen 
from  the  diagram  (fig.  4.  where  N  =  1 )  to  arise 
from  an  “identity  crisis"  in  certain  higher-order 
terms  —  i.e..  “off-circle"  terms  will  coincide  with 
“on-circle"  terms.  The  “off-circle"  coefficient 
blows  up  as  it  approaches  the  “circle." 


Also,  in  the  more  general  anisotropic  case,  even 
the  zero-order  surface  will  no  longer  be  a  circle,  but 
a  more  general  second-,  fourth-,  or  sixth-order 
surface. 


BRAGG  RESONANCE  (N=1) 


Figure  4.  An  example  of  Bragg  resonance,  which 
occurs  for  A  \sin  0  —  N  2  ithis arises analyticallv 
from  a  singularity  in  the  recursion  relation,  geo¬ 
metrically  from  the  “double  labelling"  of  an  on- 
surface  latlire  point. 

4.  COMPUTATION  AL  RESULTS 

It  is  illuminating  to  compare  results  obtained 
from  the  above  formalism  with  those  obtained 
from  earlier  efforts.  Nomoto  has  done  calculations 
for  an  isotropic  acoustically  modulated  medium 
assuming  normal  incidence.  In  this  case  (a  =  0). 
the  /-dependence  of  the  internal  fields  reduces  to 
the  tabulated  periodic  Mathieu  functions  ce„.  dis¬ 
posing  of  the  need  for  seri  ■■.<  expansions  and  so  also 
of  the  requirement  that  A  he  small.  (Nomoto 
neglects  reflections  —  that  is.  he  does  not  use  the 
boundary  conditions  —  but.  under  the  assumption 
of  a  not-too-abrupt  f-transition  at  the  interfaces, 
this  is  an  acceptable  approximation.) 

Comparisons  between  the  present  author’s  and 
Nomoto’s  calculations  are  shown  in  figure  3.  Here 
are  compared  the  intensities  of  the  undeflected 


f»  Nomofo.  Ihfbtution  of  light  hy  ullro*>*unti  f  itr’iu  ui  tht  Hnll»mn 
Ihrory,  Jit/i  f  4/i pi  l*hs*.  I U  M/  ft.’.* 
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Figure  5.  C.omparison  of  Nomoto  solution  and 
seventh-order  plane-wave  expansion;  amplitude  of 
undeflected  transmitted  beam  is  shown  as  a  func¬ 
tion  of  slab  thickness. 

^T„  =  £  Tos  exp(ip1,x)exp(iy„z)As  j 

transmitted  wave,  as  functions  of  thickness  L, 
using  the  two  methods.  For  the  system  considered. 
A  =  1 2.  The  applicability  of  the  expansion  calcu¬ 
lation.  which  was  performed  to  seventh  order,  is  in 
some  douht  for  such  a  large  A.  Also.  w/c  =  10. 
q  =  0.03.  f  =  l ,  and  L  varies  from  0  to  20  cm. 
The  value  e  =  1  was  chosen  to  minimize  the  error 
arising  from  Nomoto's  no-reflection  assumption. 
Kven  for  A  =  1 2.  the  expansion  technique  works 
well  up  to  a  thickness  of  about  10  cm  (about  500 
acoustical  wavelengths).  It  has  the  advantage  of 


applying  also  for  oblique  incidence  and  for  abrupt 
changes  in  refractive  index,  and  of  being  generaliz- 
able  to  include  anisotropic  media. 

Figure  6  shows,  for  a  different  set  of  parameter 
values,  the  amplitude  of  the  T0,  term  as  a  function 
of  incident  angle.  The  T0J  term  can  be  visualized  as 
having  been  doubly  diffracted  back  to  its  original 
direction.  Here  f  'c  =  10,  q  =  0.5.  e  =  5.  and 
L  =  1  cm.  The  Bragg  resonance  at  o  = 

(y0  =  0.25)  can  be  clearly  seen. 

The  parameter  values  in  figure  6  (e.  ca.  q.  L) 
quite  closely  resemble  those  of  the  system  that  this 
research  is  designed  to  model:  modulation  of  an 
optical  beam  by  microwave  acoustical  signals  in  a 
1-cm-thick  sample  of  lithium  niohate. 


Figure  6.  Second-order  term  of  undeflected  beam 
as  function  of  angle,  illustrating  Bragg  resonance. 
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NOTATION 


free-spaee  velocity  of  light 
optical  frequency 

dielectric  constant  of  unmodulated  material 
modulation  amplitude 

A  A2 

o  — —  .  expansion  parameter 

X 

free-spaee  optical  wavelength 
acoustic  wavelength 

'2n 

~A 


L 

y 


Rn 

Tn 

k 


slab  thickness 

z-component  of  incident  light  wave  vector 

y 

q 

term  in  expansion  of  internal  forward¬ 
travelling  wave 

term  in  expansion  of  internal  backward¬ 
travelling  wave 

term  in  expansion  of  reflected  wave 
term  in  expansion  of  transmitted  wave 
x-eomponent  of  wave  vector  in  medium 


APPENDIX  A 


APPENDIX  A.  -  THE  ACOI  STO-OPTIC  INTERACTION  FOR  BRAGG  ANGLES 


The  acousto-optic  diffraction  process  can  be 
described  by  a  perturbation  procedure  (see  main 
body  of  paper).  A  recursion  relation,  applied  in 
conjunction  with  certain  boundary  conditions, 
generates  plane-wave  expansion  terms  from  lower- 
order  terms.  The  expansion  parameter  is 


where  S  is  the  refractive  index  modulation  ampli¬ 
tude,  A  is  the  acoustic  wavelength,  and  X  is  the 
optical  wavelength.  The  recursion  relation  is 

(n  -  m)  (n  +  m  +  2«)Emns  = 

'4  (Em.n  —  l.s—  1  "E  Em  n  +  l.s— l) 

.2  s 

~~  j  2*  EmtEm  n  s_  (  (A-l) 

4ir  t  =  l 


system  unaltered.  (The  acoustical  z-period  is  an 
integral  multiple  of  the  optical  z-period.)  So  our 
solution  should  be  periodic  in  z.  After  separation  of 
variables,  the  z-part  of  the  solution  satisfies 
Mathieu's  equation;  thus,  in  the  Bragg  case,  a 
linear  combination  of  the  cen  and  sen  is  needed. 
Since  we  want  a  solution  whose  zero-order  term  is 
exp(ik„x)exp(iNqz/2),  our  solution  is  clearly 

exp(ikox)  ce\  ^-y-  ,-2a|  + 

i  exp(ik„x)  se\  ^-y-  ,-2A^  ,  (A-5) 

where 

to-  -  c  -  4- '  -  . 

C  4 


Here,  n,  m,  and  s  label  the  plane-wave  terms  which 
combine  to  give  the  full  field  pattern.  Other  vari¬ 
ables  appearing  are 

a  —  Ys‘n  &  »  (A-2) 

A 

6  —  angle  of  incidence  , 


{k-y  -  ku-=  -4- 


UJ  q2b\(-2A) 

c2  £  4 


2ttc 

~ 


2jt 


km=  L  KmlAl  , 

t=  0 

E  =  £  Emn8exp[i(kmx  +  ynz)]As  .(A-4) 

mn8 

In  the  case  of  Bragg  resonance  (a  =  N/2,  N  an 
integer),  the  method  fails,  since  the  left-hand  side 
of  (A- 1 )  vanishes.  A  consideration  of  the  solution 
for  such  cases,  and  how  to  obtain  and  represent 
them,  is  given  here. 

Physically,  the  Bragg  condition  implies  perfect 
periodicity  of  the  system;  displacement  by  A  (the 
acoustical  wavelength)  in  the  z-direction  leaves  the 


and  ajs;  and  b>,  are  the  Mathieu  eigenvalues  in 
McLachlan’s  notation.'  and  are  functions  of  A. 

This  solution  can  be  obtained  by  a  modification 
of  the  expansion  method  presented  in  the  body  of 
this  paper.  The  modification  arises  from  the  degen¬ 
eracy  (K-jv;,,  =  K,,„)  which  results  in  the  Bragg 
case.  For  definiteness,  we  consider  here  the  cast 
N  =  2  —  that  is.  y  =  q  (a  =  1).  We  use  the 
usual  recursion  relation  (A— 1),  but  we  start  with 
nonzero  K(mk,  and  E0_2I1.  Then 

Pn  £  q  “  E— 


'.V  1 1.  \h Lachlan,  Ihrory  and  -ipfilnalion  of  Malhicu  huncltons.  Oxford 
Prrss.  Oxford  (IW?) 
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For  E011  . 


The  equations  for  E^  and  E0_20  yield 


(1)(1  +  2a)E011  =  -  Eoo0 


Foil  r  Fooo 


E0-„:  (—1)1—1  +  2a)E„_„  =  \  (e^  +  E0_2o) 

Eo-ll=  2  (F°oo  '*'F0_20^ 


Eo-3i:  (  3)(  3  +  2nr)E0_j,  F0_2(l 


Fo-Jl—  ,  F0_2o 
6 


E022*  (2)(2  +  2a)E0j2  E01, 


Eo22  j6  Eon  %  Eooo 


£002:  (0)(0  +  2a)E002  =  0  = 


9  (E0ii  +E0_,|)  EojEoo,, 

q 


[M- 

1 

Eooo 

0 

q  6 

4 

1 

Eo2  1 

Eq—  20 

0 

4 

q  6j 

Thus. 


E02  V  ,  1  5 


Eo2  Eooo  2  ^  2  Eooo  2 


E0_22;  (— 2)(— 2  +  2a)E„_22  =  0  = 

~2~  (E0-ii  "hE0_3,^  —  —  K02  E0- 

K  F  _  q1  /_  1 K  _IE  ) 

^2^0-20—  2  \  2  1,000  3  c‘°~20l 

E0_42:  (-4)(-4  +  2a)E0_42=  ^E0_„ 

t’  j_  L'  _  1  L' 

E0_42=  —  Eo-.ll  —  %  E0—  20 


q’/+3  q2  144 


1  2  0  2 
E„2  =  -12  q  or  -  —  q 


Thus. 


k„  =  koo  + 


1  q‘ 


24  K 


A  +  .  .  . 


k  =  k - —  —5—  A1  + 

0  00  21  L 


Call  the  first  solution  k0  and  the  second  k„ . 
Substitution  shows  that  the  ko  solution  has  the 
property  E0_20  =  —  E000  .  and  the  k„  solution  has 
the  property  EQ_20  =  +  E000  -  We  thus  have,  to 
second  order,  the  two  solutions. 


E  =  —  E000  [exp(ik0„x)exp(iY0z) 


exp(ik00x)  exp  My, ,z)]  exp  (i 


(j  w™ ') 


2  lk00o 

=  i  E000  exp(ik0x)sin  qz  .  (A-ty 
E+  =  —  Eooo  [exp(ik00x)exp(iy0z)  + 

/  i  1  (  5A!q!  \ 

exp(ik„ox)  exp(-iy„z)Jexp  l-i  x  I 

=  Eooo  exp(ikox)cos  qz  .  ( A-7) 


14 


APPENDIX  A 


If  carried  out  to  higher  order,  these  solutions  turn 
out  to  be 

E_  =  iEooo  exp(ik^x)se2  ^-y-  ,-2a|  ,  (A-8) 
E+  =  Em,,,  exp(ik^x)  ce2  ^y  .  -2a|  . 


(A-9) 


As  remarked  earlier,  the  solution  we  want  is 

E+  +  E~  =  Egoo  ^ce2  ,-2a| 

exp(ikgx)  4-  i  se2  ^y  ,-2A^  exp(ikgx)J 

(A- 1 0) 

It  is  of  interest  now  to  investigate  the  continuity- 
properties  of  this  solution:  to  consider  the  case 
ar=l+5as5  —  0.  For  this  case,  the  reader  can 
verify  the  following: 

Kg,,  =  p0e  -  a  q  -  p0e  -  q  -  2q  0 

K_20  =  po£  -  (a  -  2)2q2  -  pot  -  q2  +  2q25  , 

Eooo  =  1  < 


F  =  — 
C-oi.  6  • 


F  =  — — 
r-o-n  2 


F  =  — 
r'°”  96 


Eo-22 


1 


0-22  7^5  * 

F  - - L 

r._,_,2  l6g 


The  E_2_22  term  represents  the  addition  of  a 
homogeneous  solution  which  can  be  considered  to 
arise  from  the  boundary  condition.  Thus,  using 

k|>a)  =  y  (k„  +  kg)  - 

1  q2  2 

k°°  12  A  ’ 
the  second-order  solution  is  given  by 

exp(ikj,a*x)  •  exp(iqz)  +  A  exp(2iqz)  —  yj 
+  ^6  exP(3iclz)  +  [65  ‘*xp(— iqz) 

(£®) 

—  exp^ikia)x^|exp(iqz)  +  A  |jr  exp(2iqz)  -  ^ 

2fl  I  iq2X 

+  A  <^exp(3iqz)-~yy-  exp(-iqz) 


(A- 11) 

This  expression  also  represents  our  Mathieu 
function  solution  to  second  order,  where  exp(ikj,a,x) 
has  been  factored  out.  The  x-dependent  coeffi¬ 
cient  arises  from  the  series  expansion  of  the 
exponentials: 


xp(ik^x) 


■f 


iexp(ik0x)  se2^y  ._2A j  = 

exp(ikj)a,x)  (l  - 
^-yjpxpl'qz)  +  exp(-iqz)j 


+ 
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A  exp(2iqz)  -  -j-  +  ^exp(-2iqz)J 

A  exP(3iqz>  +  f^2  exp(_3iqzi|  j 

exp(iki)alx)  (l  +  -^L) 


+ 


w- 


12 


exp(iqz)  —  exp( 
exp(2iqz) 


-iqz)j  + 

~  J2  exp(-2iqz)J 


+ 


'•[■IS. 


exp(3iqz) 


—  exp(-3iqz) 


•I 


The  reader  can  easily  see  that  this  matches  up  with 
the  limiting  form  of  equation  (A- 11),  and  the 
continuity  of  the  solution  at  Bragg  angles  is 
verified.  Thus,  the  expansion  procedure  described 
in  the  main  body  of  this  paper  yields  solutions  for 
the  Bragg  case  in  the  limit  a  =  N/2  +  8,  8  —  0. 
The  fact  that  the  coefficients  themselves  are  not 
continuous  arises  from  the  inherent  nonunique¬ 
ness  of  expansion  (2). 
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